Characterization and mitigation of artifacts derived from NGS library preparation due to structure-specific sequences in the human genome.

BACKGROUND: Hybridization capture-based targeted next generation sequencing (NGS) is gaining importance in routine cancer clinical practice. DNA library preparation is a fundamental step to produce high-quality sequencing data. Numerous unexpected, low variant allele frequency calls were observed in libraries using sonication fragmentation and enzymatic fragmentation. In this study, we investigated the characteristics of the artifact reads induced by sonication and enzymatic fragmentation. We also developed a bioinformatic algorithm to filter these sequencing errors. RESULTS: We used pairwise comparisons of somatic single nucleotide variants (SNVs) and insertions and deletions (indels) of the same tumor DNA samples prepared using both ultrasonic and enzymatic fragmentation protocols. Our analysis revealed that the number of artifact variants was significantly greater in the samples generated using enzymatic fragmentation than using sonication. Most of the artifacts derived from the sonication-treated libraries were chimeric artifact reads containing both cis- and trans-inverted repeat sequences of the genomic DNA. In contrast, chimeric artifact reads of endonuclease-treated libraries contained palindromic sequences with mismatched bases. Based on these distinctive features, we proposed a mechanistic hypothesis model, PDSM (pairing of partial single strands derived from a similar molecule), by which these sequencing errors derive from ultrasonication and enzymatic fragmentation library preparation. We developed a bioinformatic algorithm to generate a custom mutation "blacklist" in the BED region to reduce errors in downstream analyses. CONCLUSIONS: We first proposed a mechanistic hypothesis model (PDSM) of sequencing errors caused by specific structures of inverted repeat sequences and palindromic sequences in the natural genome. This new hypothesis predicts the existence of chimeric reads that could not be explained by previous models, and provides a new direction for further improving NGS analysis accuracy. A bioinformatic algorithm, ArtifactsFinder, was developed and used to reduce the sequencing errors in libraries produced using sonication and enzymatic fragmentation.

RNA-seq combined network pharmacology reveals that Fu-Gan-Wan (FGW) inhibits liver fibrosis via NF-κB/CCL2/CCR2 and lipid peroxidation via Nrf2/HMOX1 signaling pathway.

ETHNOPHARMACOLOGICAL RELEVANCE: Liver fibrosis is a serious complication of liver disease characterized by excessive collagen deposition, without effective therapeutic agents in the clinic. Fu-Gan-Wan (FGW) is an empirical formula used for the clinical treatment of hepatitis and cirrhosis. It has been shown to reverse experimental liver fibrosis. However, its corresponding mechanisms remain unclear. AIM OF THE REVIEW: This study aimed to elucidate the key pathways and target genes of FGW in attenuating liver fibrosis. MATERIALS AND METHODS: The therapeutic effects of different doses of FGW on liver fibrosis were investigated using a 2 mL/kg 15% CCl4-induced mouse model. Then, RNA-seq combined with network pharmacology was used to analyze the key biological processes and signaling pathways underlying the anti-liver fibrosis exertion of FGW. These findings were validated in a TGF-ß1-induced model of activation and proliferation of mouse hepatic stellate cell line JS-1. Finally, the key signaling pathways and molecular targets were validated using animal tissues, and the effect of FGW on tissue lipid peroxidation was additionally observed. RESULTS: We found that 19.5 g/kg FGW significantly down-regulated CCl4-induced elevation of hepatic ALT and AST, decreased collagen deposition, and inhibited the expression of pro-fibrotic factors α-SMA, COL1α1, CTGF, TIMP-1, as well as pro-inflammatory factor TGF-ß1. Additionally, FGW at doses of 62.5, 125, and 250 µg/mL dose-dependently blocked JS-1 proliferation, migration, and activation. Furthermore, RNA-seq identified the NF-κB signaling pathway as a key target molecular pathway for FGW against liver fibrosis, and network pharmacology combined with RNA-seq focused on 11 key genes. Significant changes were identified in CCL2 and HMOX1 by tissue RT-PCR, Western blot, and immunohistochemistry. We further demonstrated that FGW significantly attenuated CCl4-induced increases in p-p65, CCL2, CCR2, and HMOX1, while significantly elevating Nrf2. Finally, FGW significantly suppressed the accumulation of lipid peroxidation products MDA and 4-HNE and reconfigured the oxidation-reduction balance, including promoting the increase of antioxidants GPx, GSH, and SOD, and the decrease of peroxidation products ROS and GSSG. CONCLUSIONS: This study demonstrated that FGW exhibits potential in mitigating CCl4-induced hepatic fibrosis, lipid peroxidation, and iron metabolism disorders in mice. This effect may be mediated through the NF-κB/CCL2/CCR2 and Nrf2/HMOX1 pathways.

Astragaloside IV targeting autophagy of T cells improves inflammation of asthma.

Astragaloside IV (AST) has been confirmed to have antiasthmatic effects. However, the underline mechanism is unclear. The study aimed to explore the treatment mechanism of AST based on autophagy of memory T cells. AST treatment significantly decreased the number of T effector cells in asthma mice blood and the nude mice that received AST-treated TCMs had relieved inflammation compared with the untreated group; meanwhile, we found that AST significantly decreased the autophagy level and inhibited OX40/OX40L signal pathway of lymphocytes. The results highlighted that AST regulated autophagy to inhibit differentiation of effector T-cell phenotype.

Atp6v1h Deficiency Blocks Bone Loss in Simulated Microgravity Mice through the Fos-Jun-Src-Integrin Pathway.

The microgravity conditions in outer space are widely acknowledged to induce significant bone loss. Recent studies have implicated the close relationship between Atp6v1h gene and bone loss. Despite this, the role of Atp6v1h in bone remodeling and its molecular mechanisms in microgravity have not been fully elucidated. To address this, we used a mouse tail suspension model to simulate microgravity. We categorized both wild-type and Atp6v1h knockout (Atp6v1h+/-) mice into two groups: regular feeding and tail-suspension feeding, ensuring uniform feeding conditions across all cohorts. Analysis via micro-CT scanning, hematoxylin-eosin staining, and tartrate-resistant acid phosphatase assays indicated that wild-type mice underwent bone loss under simulated microgravity. Atp6v1h+/- mice exhibited bone loss due to Atp6v1h deficiency but did not present aggravated bone loss under the same simulated microgravity. Transcriptomic sequencing revealed the upregulation of genes, such as Fos, Src, Jun, and various integrin subunits in the context of simulated microgravity and Atp6v1h knockout. Real-time quantitative polymerase chain reaction (RT-qPCR) further validated the modulation of downstream osteoclast-related genes in response to interactions with ATP6V1H overexpression cell lines. Co-immunoprecipitation indicated potential interactions between ATP6V1H and integrin beta 1, beta 3, beta 5, alpha 2b, and alpha 5. Our results indicate that Atp6v1h level influences bone loss in simulated microgravity by modulating the Fos-Jun-Src-Integrin pathway, which, in turn, affects osteoclast activity and bone resorption, with implications for osteoporosis. Therefore, modulating Atp6v1h expression could mitigate bone loss in microgravity conditions. This study elucidates the molecular mechanism of Atp6v1h's role in osteoporosis and positions it as a potential therapeutic target against environmental bone loss. These findings open new possibilities for the treatment of multifactorial osteoporosis.

Potential Theranostic Roles of SLC4 Molecules in Human Diseases.

The solute carrier family 4 (SLC4) is an important protein responsible for the transport of various ions across the cell membrane and mediating diverse physiological functions, such as the ion transporting function, protein-to-protein interactions, and molecular transduction. The deficiencies in SLC4 molecules may cause multisystem disease involving, particularly, the respiratory system, digestive, urinary, endocrine, hematopoietic, and central nervous systems. Currently, there are no effective strategies to treat these diseases. SLC4 proteins are also found to contribute to tumorigenesis and development, and some of them are regarded as therapeutic targets in quite a few clinical trials. This indicates that SLC4 proteins have potential clinical prospects. In view of their functional characteristics, there is a critical need to review the specific functions of bicarbonate transporters, their related diseases, and the involved pathological mechanisms. We summarize the diseases caused by the mutations in SLC4 family genes and briefly introduce the clinical manifestations of these diseases as well as the current treatment strategies. Additionally, we illustrate their roles in terms of the physiology and pathogenesis that has been currently researched, which might be the future therapeutic and diagnostic targets of diseases and a new direction for drug research and development.

A unified DNA- and RNA-based NGS strategy for the analysis of multiple types of variants at the dual nucleic acid level in solid tumors.

BACKGROUND: Targeted next-generation sequencing (NGS) is a powerful and suitable approach to comprehensively identify multiple types of variants in tumors. RNA-based NGS is increasingly playing an important role in precision oncology. Both parallel and sequential DNA- and RNA-based approaches are expensive, burdensome, and have long turnaround times, which can be impractical in clinical practice. A streamlined, unified DNA- and RNA-based NGS approach is urgently needed in clinical practice. METHODS: A DNA/RNA co-hybrid capture sequencing (DRCC-Seq) approach was designed to capture pre-capture DNA and RNA libraries in a single tube and convert them into one NGS library. The performance of the DRCC-Seq approach was evaluated by a panel of reference standards and clinical samples. RESULTS: The average depth, DNA data ratio, capture ratio, and target coverage 250 (×) of the DNA panel data had a negative correlation with an increase in the proportion of RNA probes. The SNVs, indels, fusions, and MSI status were not affected by the proportion of RNA probes, but the copy numbers of the target genes were higher than expected in the standard materials, and many unexpected gene amplifications were found using D:R (1:2) and D:R (1:4) probe panels. The optimal ratio of DNA and RNA probes in the combined probe panel was 1:1 using the DRCC-Seq approach. The DRCC-Seq approach was feasible and reliable for detecting multiple types of variants in reference standards and real-world clinical samples. CONCLUSIONS: The DRCC-Seq approach is more cost-effective, with a shorter turnaround time and lower labor requirements than either parallel or sequential targeted DNA NGS and RNA NGS. It is feasible to identify multiple genetic variations at the DNA and RNA levels simultaneously in clinical practice.

Genotype-phenotype analysis of selective failure of tooth eruption-A systematic review.

Tooth eruption is an important and unique biological process during craniofacial development. Both the genetic and environmental factors can interfere with this process. Here we aimed to find the failure pattern of tooth eruption among five genetic diseases. Both systematic review and meta-analysis were used to identify the genotype-phenotype associations of unerupted teeth. The meta-analysis was based on the characteristics of abnormal tooth eruption in 223 patients with the mutations in PTH1R, RUNX2, COL1A1/2, CLCN7, and FAM20A respectively. We found all the patients presented selective failure of tooth eruption (SFTE). Primary failure of eruption patients with PTH1R mutations showed primary or isolated SFTE1 in the first and second molars (59.3% and 52% respectively). RUNX2 related cleidocranial dysplasia usually had SFTE2 in canines and premolars, while COL1A1/2 related osteogenesis imperfecta mostly caused SFTE3 in the maxillary second molars (22.9%). In CLCN7 related osteopetrosis, the second molars and mandibular first molars were the most affected. While FAM20A related enamel renal syndrome most caused SFTE5 in the second molars (86.2%) and maxillary canines. In conclusion, the SFTE was the common characteristics of most genetic diseases with abnormal isolated or syndromic tooth eruption. The selective pattern of unerupted teeth was gene-dependent. Here we recommend SFTE to classify those genetic unerupted teeth and guide for precise molecular diagnosis and treatment.

Molecular Mechanisms of Craniofacial and Dental Abnormalities in Osteopetrosis.

Osteopetrosis is a group of genetic bone disorders characterized by increased bone density and defective bone resorption. Osteopetrosis presents a series of clinical manifestations, including craniofacial deformities and dental problems. However, few previous reports have focused on the features of craniofacial and dental problems in osteopetrosis. In this review, we go through the clinical features, types, and related pathogenic genes of osteopetrosis. Then we summarize and describe the characteristics of craniofacial and dental abnormalities in osteopetrosis that have been published in PubMed from 1965 to the present. We found that all 13 types of osteopetrosis have craniomaxillofacial and dental phenotypes. The main pathogenic genes, such as chloride channel 7 gene (CLCN7), T cell immune regulator 1 (TCIRG1), osteopetrosis-associated transmembrane protein 1 (OSTM1), pleckstrin homology domain-containing protein family member 1 (PLEKHM1), and carbonic anhydrase II (CA2), and their molecular mechanisms involved in craniofacial and dental phenotypes, are discussed. We conclude that the telltale craniofacial and dental abnormalities are important for dentists and other clinicians in the diagnosis of osteopetrosis and other genetic bone diseases.

Molecular-based phenotype variations in amelogenesis imperfecta.

Amelogenesis imperfecta (AI) is one of the typical dental genetic diseases in human. It can occur isolatedly or as part of a syndrome. Previous reports have mainly clarified the types and mechanisms of nonsyndromic AI. This review aimed to compare the phenotypic differences among the hereditary enamel defects with or without syndromes and their underlying pathogenic genes. We searched the articles in PubMed with different strategies or keywords including but not limited to amelogenesis imperfecta, enamel defects, hypoplastic/hypomaturation/hypocalcified, syndrome, or specific syndrome name. The articles with detailed clinical information about the enamel and other phenotypes and clear genetic background were used for the analysis. We totally summarized and compared enamel phenotypes of 18 nonsyndromic AI with 17 causative genes and 19 syndromic AI with 26 causative genes. According to the clinical features, radiographic or ultrastructural changes in enamel, the enamel defects were basically divided into hypoplastic and hypomineralized (hypomaturated and hypocalcified) and presented a higher heterogeneity which were closely related to the involved pathogenic genes, types of mutation, hereditary pattern, X chromosome inactivation, incomplete penetrance, and other mechanisms.The gene-specific enamel phenotypes could be an important indicator for diagnosing nonsyndromic and syndromic AI.

Hereditary dentin defects with systemic diseases.

OBJECTIVE: This review aimed to summarize recent progress on syndromic dentin defects, promoting a better understanding of systemic diseases with dentin malformations, the molecules involved, and related mechanisms. SUBJECTS AND METHODS: References on genetic diseases with dentin malformations were obtained from various sources, including PubMed, OMIM, NCBI, and other websites. The clinical phenotypes and genetic backgrounds of these diseases were then summarized, analyzed, and compared. RESULTS: Over 10 systemic diseases, including osteogenesis imperfecta, hypophosphatemic rickets, vitamin D-dependent rickets, familial tumoral calcinosis, Ehlers-Danlos syndrome, Schimke immuno-osseous dysplasia, hypophosphatasia, Elsahy-Waters syndrome, Singleton-Merten syndrome, odontochondrodysplasia, and microcephalic osteodysplastic primordial dwarfism type II were examined. Most of these are bone disorders, and their pathogenic genes may regulate both dentin and bone development, involving extracellular matrix, cell differentiation, and metabolism of calcium, phosphorus, and vitamin D. The phenotypes of these syndromic dentin defects various with the involved genes, part of them are similar to dentinogenesis imperfecta or dentin dysplasia, while others only present one or two types of dentin abnormalities such as discoloration, irregular enlarged or obliterated pulp and canal, or root malformation. CONCLUSION: Some specific dentin defects associated with systemic diseases may serve as important phenotypes for dentists to diagnose. Furthermore, mechanistic studies on syndromic dentin defects may provide valuable insights into isolated dentin defects and general dentin development or mineralization.

[Research progress on structure, function and disease correlation of solute carrier family 4].

The solute carrier family 4 (SLC4) includes 10 members (SLC4A1-5, SLC4A7-11), which are expressed in multiple tissues in the human body. The SLC4 family members differ in their substrate dependence, charge transport stoichiometry and tissue expression. Their common function is responsible for the transmembrane exchange of multiple ions, which is involved in many important physiological processes, such as erythrocyte CO2 transport and the regulation of cell volume and intracellular pH. In recent years, many studies have focused on the role of SLC4 family members in the occurrence of human diseases. When SLC4 family members have gene mutations, a series of functional disorders will occur in the body, leading to the occurrence of some diseases. This review summarizes the recent progress about the structures, functions and disease correlation of SLC4 members, in order to provide clues for the prevention and treatment of related human diseases.

Baicalin Inhibits Airway Smooth Muscle Cells Proliferation through the RAS Signaling Pathway in Murine Asthmatic Airway Remodeling Model.

Background: Studies that looked at asthma airway remodeling pathogenesis and prevention have led to the discovery of the rat sarcoma viral oncogene (RAS) signaling pathway as a key mechanism that controls airway smooth muscle cell (ASMC) proliferation. Baicalin has great anti-inflammatory, proliferation-inhibited, and respiratory disease-relieving properties. However, the inhibitory effects and mechanisms of baicalin on ASMC-mediated airway remodeling in mice are still poorly understood. Methods: After establishing the asthmatic mice model by ovalbumin (OVA) and interfering with baicalin, airway remodeling characteristics such as airway resistance, mRNA, and protein expression levels of remodeling-related cytokines were measured by histopathological assessment, quantitative real-time polymerase chain reaction (qPCR), enzyme-linked immunosorbent assay (ELISA), and western blot. Further efforts on detailed mechanisms were used antibody arrays to compare the expression and activation of proteins involved in the RAS signaling pathway. In addition, validation experiments were performed in ASMC proliferation model and low-expression cells of the target gene by using shRNA. Results: In OVA-induced asthmatic mice model, baicalin significantly reduced the infiltration of inflammatory cells in lung tissue, attenuated airway resistance, and decreased mRNA and protein expression levels of remodeling-related cytokines such as interleukin-13 (IL-13), vascular endothelial growth factor (VEGF), transforming growth factor-beta 1 (TGF-ß1), matrix metallopeptidase 9 (MMP9), and tissue inhibitor of metalloproteinase 1 (TIMP1). The results of antibody arrays involved in RAS signaling pathway revealed that OVA and baicalin administration altered the activation of protein kinase C alpha type (PKC-α), A-rapidly accelerated fibrosarcoma (A-RAF), mitogen-activated protein kinase 2 (MEK2), extracellular regulated MAP kinase (ERK), MAPK interacting serine/threonine kinase 1 (MNK1), and ETS transcription factor 1 (ELK1). The above results were further verified in the ASMC proliferation model. A-RAF silencing (shA-RAF) could promote ASMC proliferation and downregulate p-MEK2, p-ERK, p-MNK1, and p-ELK1 expression. Conclusion: The effects of baicalin against airway remodeling and ASMC proliferation might partially be achieved by suppressing the RAS signaling pathway. Baicalin may be a new therapeutic option for managing airway remodeling in asthma patients.

Case report: identification of EGFR R776H and FANCE R381H germline mutations in a patient with multiple pulmonary nodules.

Pulmonary nodules evaluation is clinically crucial because they may be the early predictors of lung cancer. Except for CT screening and serum tumor biomarkers testing, genetic alteration analysis by next-generation sequencing (NGS) technology can also help to find cancer earlier. In this study, we report a case of multiple pulmonary nodules patient with EGFR R776H and FANCE R381H germline mutations. Her father, paternal aunt, and elder uncle harbored either one or both two mutations and were found with multiple pulmonary ground-glass or sub-solid nodules. Her 7-year-old daughter also inherited the same two mutations.

MSX1 involved selective tooth agenesis and abnormal labial frenum, pedigree, and retrospective study.

OBJECTIVE: Muscle segment homeobox gene 1 (MSX1) is widely expressed in craniofacial development and tooth formation. The aim of this study was to report a novel MSX1 mutation in a Chinese family with selective tooth agenesis and abnormal median maxillary labial frenum (MMLF). MATERIALS AND METHODS: Mutation analysis was carried out by whole exome sequencing. The pMD18-T vector was used to verify the mutations. PubMed and Human Gene Mutation Database were searched to analyze the relationship between the mutations in MSX1 and related phenotypes. RESULTS: A novel heterozygous mutation (c.75delG) in MSX1 was detected in the proband and her mother. They presented as oligodontia and lower attached hypertrophy median maxillary labial frenum. 60 MSX1 mutations from 39 reports did not declare malformed MMLF except our cases. Meanwhile, we found that the types and sites of MSX1 mutations may affect the selectivity of tooth agenesis and orofacial cleft. CONCLUSION: This study suggests malformed MMLF as a new phenotype of MSX1 mutation and a specific relationship between MSX1 genotype and phenotype.

Exosomal microRNA panel as a diagnostic biomarker in patients with hepatocellular carcinoma.

Background: Diagnostic tools for hepatocellular carcinoma (HCC) are critical for patient treatment and prognosis. Thus, this study explored the diagnostic value of the exosomal microRNA panel for HCC. Methods: Expression profiles of microRNAs in exosomes and plasma of HCC and control groups were assessed using microRNA microarray analysis. Reverse transcription-quantitative PCR was applied to evaluate the expression of candidate microRNAs in blood samples from 50 HCC patients, 50 hepatic cirrhosis patients, and 50 healthy subjects. The area calculated the diagnostic accuracy of the microRNAs and microRNA panel under the receiver operating characteristic curve (AUC). Results: MicroRNA microarray analysis revealed that there were more differentially expressed microRNAs in the exosome HCC group than plasma HCC group. Among the 43 differentially expressed microRNAs contained in both exosomes and plasma, we finally decided to testify the expression and diagnostic significance of microRNA-26a, microRNA-29c, and microRNA-199a. The results indicated that expression of the microRNA-26a, microRNA-29c, and microRNA-199a in both exosomes and plasma was significantly lower in HCC patients compared with hepatic cirrhosis and healthy group. Interestingly, exosomal microRNAs were substantially more accurate in diagnosing HCC than microRNAs and alpha-fetoprotein in plasma. Moreover, the exosomal microRNA panel containing microRNA-26a, microRNA-29c, and microRNA-199a showed high accuracy in discriminating HCC from healthy (AUC = 0.994; sensitivity 100%; specificity 96%) and hepatic cirrhosis group (AUC = 0.965; sensitivity 92%; specificity 90%). Conclusion: This study revealed that the exosomal microRNA panel has high accuracy in diagnosing HCC and has important clinical significance.

OncoPubMiner: a platform for mining oncology publications.

Updated and expert-quality knowledge bases are fundamental to biomedical research. A knowledge base established with human participation and subject to multiple inspections is needed to support clinical decision making, especially in the growing field of precision oncology. The number of original publications in this field has risen dramatically with the advances in technology and the evolution of in-depth research. Consequently, the issue of how to gather and mine these articles accurately and efficiently now requires close consideration. In this study, we present OncoPubMiner (https://oncopubminer.chosenmedinfo.com), a free and powerful system that combines text mining, data structure customisation, publication search with online reading and project-centred and team-based data collection to form a one-stop 'keyword in-knowledge out' oncology publication mining platform. The platform was constructed by integrating all open-access abstracts from PubMed and full-text articles from PubMed Central, and it is updated daily. OncoPubMiner makes obtaining precision oncology knowledge from scientific articles straightforward and will assist researchers in efficiently developing structured knowledge base systems and bring us closer to achieving precision oncology goals.

K205R specific nanobody-horseradish peroxidase fusions as reagents of competitive ELISA to detect African swine fever virus serum antibodies.

BACKGROUND: African swine fever virus (ASFV) is a highly contagious hemorrhagic disease and often lethal, which has significant economic consequences for the swine industry. Due to lacking of commercial vaccine, the prevention and control of ASF largely depend on early large-scale detection and screening. So far, the commercial ELISA kits have a long operation time and are expensive, making it difficult to achieve large-scale clinical applications. Nanobodies are single-domain antibodies produced by camelid animals, and have unique advantages such as smaller molecular weight, easy genetic engineering modification and low-costing of mass production, thus exhibiting good application prospects. RESULTS: The present study developed a new method for detection of ASFV specific antibodies using nanobody-horseradish peroxidase (Nb-HRP) fusion proteins as probe. By using camel immunization, phage library construction and phage display technology, five nanobodies against K205R protein were screened. Then, Nb-HRP fusion proteins were produced using genetic modification technology. Based on the Nb-HRP fusion protein as specific antibodies against K205R protein, a new type of cELISA was established to detect ASFV antibodies in pig serum. The cut-off value of the cELISA was 34.8%, and its sensitivity, specificity, and reproducibility were good. Furthermore, the developed cELISA exhibited 99.3% agreement rate with the commercial available ELISA kit (kappa value = 0.98). CONCLUSIONS: The developed cELISA method has the advantages of simple operation, rapid and low-costing, and can be used for monitoring of ASFV infection in pigs, thus providing a new method for the prevention and control of ASF.

Molecular Features in Lymphatic Metastases Reflect the Metastasis Mechanism of Lymph Nodes With Non-Small-Cell Lung Cancers.

Lymphatic metastasis influences clinical treatment and prognosis of patients with non-small-cell lung cancer (NSCLC). There is an urgency to understand the molecular features and mechanisms of lymph node metastasis. We analyzed the molecular features on pairs of the primary tumor and lymphatic metastasis tissue samples from 15 NSCLC patients using targeted next-generation sequencing. The potential metastasis-related genes were screened from our cohort based on cancer cell fraction. After filtering with gene functions, candidate metastasis-related events were validated in the MSK cohort with Fisher's exact test. The molecular signature and tumor mutational burden were similar in paired samples, and the average mutational concordance was 42.0% ± 28.9%. Its metastatic mechanism is potentially a linear progression based on the metastatic seeding theory. Furthermore, mutated ataxia telangiectasia mutated and Rad3-related (ATR) and tet methylcytosine dioxygenase 2 (TET2) genes were significantly enriched in lymphatic metastases (p ≤ 0.05). Alterations in these two genes could be considered metastasis-related driving events. Mutated ATR and TET2 might play an active role in the metastasis of lymph nodes with NSCLC. More case enrollment and long-term follow-up will further verify the clinical significance of these two genes.